pub fn new_birthday_probability(n: u32) -> f64 {
    if n > 365 {
        return 1.;
    }
    1. - ((365 - n + 1)..=365).map(f64::from).product::<f64>() / 365f64.powi(n as _)
}
